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Abstract. In this review, we examine the successes and weaknesses of modern 
low-mass star and brown dwarf theory from various comparisons with available 
experimental and observational constraints in different domains. (1) We first 
focus on the mechanical (equation of state) and thermal (atmosphere) properties 
and on the evolution of such cool and dense objects. We then examine the current 
shortcomings of the theory and we discuss in detail recent observational analysis 
which have suggested discrepancies between the models and the observations. (2) 
We then examine the stellar and brown dwarf initial mass function and suggest 
that a power-law above the average thermal Jeans mass (~ 1M©) rolling over a 
lognormal form below this limit adequately reproduces the observations of field 
and young cluster stellar and brown dwarf distributions. This yields a reasonably 
accurate estimate of the stellar and brown dwarf Galactic census. Finally (3) we 
examine the modern context of star and brown dwarf formation and argue that 
the combination of turbulence driven fragmentation at large scale and gravity at 
small scales provides an appealing solution for the general star and brown dwarf 
formation mechanism. It also provides a physical ground for the aforementioned 
power-law + lognormal form for the IMF, whereas a series of different power laws 
lacks such a physical motivation. Finally we argue that the deuterium-burning 
limit as the distinction between stars and planets has no physical foundation in 
this modern star formation scheme and should be abandoned. Opacity limited 
fragmentation extending below 10 jupiter-masses down to a few jupiter masses, 
due to shocks, anisotropy or magnetic fields, provides a much more robust limit, 
even though difficult to determine accurately. Therefore, the various "direct" 
detections of exoplanets claimed recently in the literature are most likely regular 
low-mass brown dwarfs and the direct detection of an extrasolar planet remains 
for now elusive. 



1. Introduction 

Ongoing observational projects have revealed hundreds of substellar objects 
(SSO), brown dwarfs (BD) and gaseous planets orbiting stars outside the so- 
lar system. These surveys deliver large observational data bases and necessitate 
the best possible theoretical foundation not only to understand the properties of 
these objects, including their formation mechanism, but also to evaluate prop- 
erly their contributions to the baryonic content of the Galaxy. Our group has 
ambitioned to derive a complete theory of the structure, evolution and spectral 
signature of low-mass, dense objects, from Sun-like to Saturn-like masses, cover- 
ing 3 orders of magnitude in mass and 9 in luminosity, bridging the gap between 
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stars and gaseous planets. We have incorporated the best possible physics aimed 
at describing the mechanical and thermal properties of these objects - equation 
of state, synthetic spectra, non-grey atmosphere models. We refer to previously 
published papers and reviews for a complete description of the physics entering 
these models. In this review we briefly examine the present status of this theory 
by confronting the theroretical predictions to various experimental and obser- 
vational results. We then explore the Galactic implications of the physics of 
low-mass objects, namely mass function and Galactic census. Finally, we exam- 
ine the formation of brown dwarfs and low-mass stars in the context of modern 
star formation theory. 



2. The physics of low- mass stars and brown dwarfs 
2.1. Interior 

Central conditions for low- mass stars (LMS), i.e. objects below one solar mass, 
and brown dwarfs, i.e. objects below the hydrogen-burning minimum mass 
(mHBMM ^ 0.07 M ), range from T c ~ 10 4 -10 7 K and p c ~ 10-10 3 g cm" 3 , 
when spaning conditions from the Sun to Jupiter- mass interiors. Under these 
conditions, the average ion electrostatic energy (Ze) 2 /a, where a = (^:^-) 1/ ' 3 
is the mean interionic distance, is several times the average kinetic energy kT, 
characterizing a strongly coupled ion plasma. The temperature is of the order 
of the electron Fermi temperature kTp, so that the electron gas is only partially 
degenerate. The temperature in the envelope is kT < 1 Ryd, so electronic and 
atomic recombinations take place in this region. Finally, the electron binding 
energy is of the order of the Fermi energy Ze 2 /ao ~ Ep so that pressure dis- 
sociation and ionization take place along the interior profile. Within the past 
recent years, several high-pressure shock wave experiments have been conducted 
in order to probe the equation of state (EOS) of deuterium, the isotope of hy- 
drogen, under conditions characteristic of the interior of these objects. Gas 
gun shoc k compression exp eriments were generally limited to pressures below 
1 Mbar ( Nellis et al. 1983 ). probing only the domain of molecular hydrogen. 



The Saumon-Chabrier-VanHorn EOS (Saumon, Chabrier & VanHorn 1995) was 
found to adequately reproduce the experimental pressure-density profile but to 
yield temperatures about 30% higher than the gas gun results, indicating in- 
suffisant H2 dissociati on. A slightly revise d version recovers well these gas gun 
experimental results ( Saumon et al. 2000t) . Moder n techniques include laser- 



driven shock - wave experiments dCollins et al. 1998T ). pulse-power compression 
experi ments dKnudson et al. 2 004 ) and conver gent spherical shock wave experi- 
ments (|Belov et al. 20021 : iBoriskov et al."2 003) and can achieve higher pressures 



than gas gun experiments, namely up to 5 Mbar on liquid D2 at high temper- 
ature, exploring for the first time the regime of pressure-dissociation and ion- 
ization. Unfortunately, these recent experiments give different results at high 
pressure and this controversy remains to be settled before robust comparison 
between experiment and theory can be made in the very domain of hydrogen 
pressure-ionization, i.e. P ~ 1-3 Mbar. For pressures above a few Mbar, hy- 
drogen and helium become fully ionized and accurate Monte Carlo simulations 
of the thermodynamic properties can be used as a benchmark for the EOS. 
These high-pressure experiments clearly open a new window in physics and as- 
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trophysics by constraining the EOS of stellar, brown dwarf and planet interiors 
in laboratory experiments. 

2.2. Atmosphere 

Atmosphere equilibrium condition dr = RdP/g, where g = Gm/R 2 ~ 10 3 — 10 5 
cm s -2 is the characteristic surface gravity of SSOs, yields P p h ~ g/R « 0.01 — 1 
bar and p p h «10~ 6 -10~ 5 g cnr 3 near the photosphere. Collision effects are 
important under such conditions and induce new sources of absorption like 
e.g. the collision-induced absorption between H2-H2 or H2-He below ~ 5000 
K ( Linskv 19691 : iBorvsow et al. 19 85). For effective temperatures characteris- 



tic of LMS and BDs (T e g<5000 K), numerous molecules form, in particular 
metal oxydes and hydrides (TiO, VO, FeH, CaH, MgH), the major absorbers 
in the optical, and CO, H2O which dominate in the infrared. The situation 
becomes even more complex for temperatures below T ~ 2500 K, i.e. for atmo- 
spheric conditions of objects ranging from the coolest stars to BDs and jovian 
planets. For T e g- < 2500 K, i.e. at the bottom of the main sequence, there is 
evidence for condensation of metals and silicates into grai ns (e.g. TiO into 
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into carbon monoxide CO, while the oxygen is found in titanium TiO and vana- 
dium VO monoxides (dominating the optical absorption) and water vapor H2O 
(shaping the infrared spectrum). This absence of refractory metal line absorp- 
tion signature in the spectrum defines the so-called "L-dwarf" domain. Below 
y 1800 K, the dominant equilibrium form of carbon is no longer CO but CH4 
( Fegl ev fc Lodders 19 96 ) , as identified f or e.g. in G1229B (|Qppenheimer et al. 1995; 



lAllard et al. 19961: I Marlev et al. 19961 ) . As confirmed by the observation of 
G1229B (|Qppenheimer et al. 19951 ). methane bands begin to appear in the in- 
frared (1.6 and 2.2 /im), while titanium dioxide and silicate clouds form at the 
expense of TiO, modifying profoundly the thermal opacity of the atmosphere. 
The presence of methane absorption in the spectrum defines the "T-dwarf ' spec- 
tral type, even though weak CH4 absorption is already present in the latest L 
dwarfs. For jovian-like atmospheres, the dominant equilibrium form of nitrogen 
N2 is NH3 (T e g <; 600 K) and below T e fj ~ 600 K (depending on gravity) water 
condenses to clouds at or above the photosphere. Figures Q displays the spectral 
evolution of M, L and T dwarfs, with identified main absorption sources. In 
spite of all these complex phenomena, BDs radiate nearly 90% of their energy 
at wavelengths longward of 1 /mi, with a peak around ~ 1-2 /mi, and infrared 
broad-band colors are prefered to optical ones for the identification of these 
objects. 

The condensates or grains affect the atmosphere in different ways. Grain 
formation depletes the corresponding gas-phase absorber in certain regions of the 
atmosphere and modifies the EOS itself and thus the atmospheric temperature 
profile, but it also strongly modifies the opacities and thus the emergent spec- 
trum. At last it produces an increase of the temperature in the uppermost layers 
of the atmosphere, where lines form, the backwarming effect, so that molecular 
bands like e.g. H2O form in hotter regions and are thus weaker. This produces 
a significant reddening of the colors in late M and L dwarfs. In the opposite, the 
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Figure 1. Spectral energy distributions of typical M, L and T dwarfs ac- 
cording to models of Allard et al. (2001). From top to botton, T c s — 2500, 
1800 and 1000 K, respectively, with surface gravities and radii from Chabrier 
ct al. (2000) for an age 5 Gyr. The most important bands of absorption are 
indicated. 



strong absorption bands of methane in the infrared (H, K and L bands) induce 
a redistribution of the emergent flux at shorter wavelengt hs, yielding bluer co l- 
ors. Our previous calculations of main sequence objects ([Baraffe at al. 1 998). 
based on the atmosphere models of Hauschildt et al. (1999) have been ex- 
tended by taking into account this grain formation process in two extreme 
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regimes. In the so-called "dusty" limit, all condensed species are included ex- 
plicitely both in th e atmosphere EOS and in the radiative transfer equation 
((Allard et al. 2001ft and provide the proper boun dary conditions for con sistent 
evolutionary calculations for such "dusty" objects ( Chabrier et al. 2000aft in the 
effective temperature range 1700 ^ T e g- ^ 2500 K. At low enough temperature, 
T e ff <^ 1700 K, we consider the other so-called "condensed" limit, where all 
the grains either form or have sunk below the photos phere, halting the red- 
ward progression of infrared colors ( Baraffe et al. 20031 1. Several objects have 
now been discovered in the in-between "early T-dwarfs" domain, with 0.5 ^ J- 
K <^ 1.5 (Kir kpatrick et al. 20001 : iBurgasser et al. 20 02). In this regime, grain 
dynamics with competing cqnvective and settlin g time scales must be considered 
( Ackerman fe Marlev 200 3 : 1 Allard et al. 2003aft . yielding complex, and possibly 
non-mono tonic color-distributions with brightening J luminosities fo r the earliest 
T dwarfs ((Tinnev et al. 20031 : lKn"app et al. 2003 : IVrba et al. 2004ft . 

The "dusty" models have been shown to reproduce the spectra and colors 
of various early L-dwarfs including t he ones discovered by the DENIS project 
( Ruiz et al. 19971 : iTinney et al. 1998ft and GD16 5B. which lies on the very edge 
of the H-burning limit (iKirknatrick et al. 1 9991 while the "condensed' mod- 
els give a go od representation of the T-dwarf spect ral energy distribution and 
photometry ((Allard et al. 199fil : iBaraffe et al. 2003ft . In the domain between 
late L-dwarfs and T-dwarfs, as mentioned earlier, models including dynamical 
processes must be considered. Atmosphere models including cloud sedimen- 
tation seem to give at least a qualitatiy e description of the L to T transition 
( Marlev et al. 20021 : lAllard et al. 2003aft . No consistent evolutionary calcula- 
tion, however, exists in this domain yet. Fortunately, this late-L early-T dwarf 
domain concerns only massive and/or young BDs, and thus a limited fraction of 
the global BD population. 

Although the general description of BD spectral evolution can be considered 
as satisfactory, there is definitely room for improvement in a few bands, where 
the observed flux is either overestimated or underestimated by the models. Such 
discrepancies are mostly tied to incomplete opacities (e.g. TiO, H2O, CH4, CaH, 
VO, FeH). An important issue in the T-dwarf domain concerns the alkali- line 

pressure broadening like e.g. the Na I (0.59 jum) and K I (0.77 ^m) resonance 

doublets ((Tsuii et al. 19991 : IBurrows et al. 2000l : iBurrows fc Volobuvev 20031 : lAllard et al. 2003bft 
which absorb the flux in the optical, yieding redder optical colors with de- 
creasing T e g, and accurate treatments of the far wings of these lines are nec- 
essary to model accurately very cool BDs. The spectral evolution of BDs is 
rendered even more complicated by the occurence of non equilibrium ther- 
mochemical reactions, as observed e.g. by the greater absorption of CO and 
thus an enhanced abundance than pr edicted by chemical equilibrium models 
( Noll et al. 19971 : l aumon et al. 2003ft . At last, there is strong observational 
suggestion that the cloud surface coverage breaks up in the L-T transition re- 
gion, as suggested by the aforementioned incr easing J flux, which c ould then 
emerge from deeper layers of the atmosphere ( Burgasser et al . 2002). Spectra 
changes in this domain, i.e. around ~ 1300 K, are thus not due to decreasing 
temperature, this latter remaining constant over the ~ L7-T4 domain, but coul d 
be explained by clearing of the cloud decks (|Burgasser et al. 2 002: Tsuii 2005). 
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3. Evolution 

3.1. Color-magnitude diagrams of disk and young cluster objects 




Figure 2. Ncar-IR K — (J — K) color-magnitude diagram. The data are 
from Leggett (1992) (dots), Dahn et al. (2000) (squares) and Reid et al. 
(1999) (filled triangles). Theoretical models correspond to various atmosphere 
models and ages: DUSTY (right) and COND (left) models. Preliminary 
models for late L early T are also shown between these two limits. Masses (in 
Mq) and T cS are indicated for the DUSTY and COND isochrones. 

As mentioned above, present evolutionary calculations have been conducted 
with three different atmosphere models: (i) for objects hot enough, i.e. massive 
or young enough, to preclude the formation of grains (T e g ^2800 K), we have 
used the grainless so -called "NextGen" atmosphere models (|Hauschildt et al. 19991 : 
iBaraffe at al. 1998ft : (ii) for objects below 2800 K, i.e. the L-dwarf regime, 
which encompasses the bottom of the main sequence, we have used atmo- 
sphere models wich include the grain opacity sour ces in the transfer equ ation 
( Allard et al. 200lT l. the so-called "dusty" models (jChabrier et al. 200 03): (hi) 
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for objects below about 1500 K, down to Jupiter-like temperatures, we have con- 
sidered cases where the condensates settle rapidly below the photosphere and 
- although modifying the atmosphere E OS - do not particip ate to the opacity, 
the aforementioned "condensed" models (|Baraffe et al. 20 03 ) . These latter cal- 
culations are motivated by the absence of grain features in the atmosphere of 
objects below T e g ~ 1000 K, i.e. Gliese 229B-like objects. 

Figure 21 displays a K vs J-K color-magnitude diagram (CMD). In terms 
of colors there is a competing effect between grain and molecular sources of 
absorption for objects at the bottom of and just below the main sequence, with 
grain opacity yielding eventually a severe redshift of the colors, a consequence 
of the aforementioned backwarming effect. For cooler objects, CH4 absorption 
in the infrared leads to a blueshift for infrared colors. The figure also displays 
preliminary calculat ions based on atmosphere models including grain dynamics 
( Allard et al. 2003al ) in the late-L early-T domain. 

Several BD surveys have been conducted in young clusters and it is im- 
portant to develop accurate non-grey models for PMS stars and young BDs. 
The more massive of such objects are too hot for grain formation to occur. Al- 
though the Baraffe et al . (1998) models giv e a consistent description of the 
CMDs of these clusters (|Luhman et al. 2003 i. the situation for young objects 
is not as satisfactory as for objects with higher gravity. As discussed in detail 
in Baraffe et al. (2002) both models and observations are hampered by numer- 
ous uncertainties and great caution must be taken when considering young age 
(<; 10 Myr) objects. Note that, when using the Luhman (1999) T e g-spectral 
type scale for young objects adjusted to yield consistent, coeval sequences for 
the quadruple system GG-tau, whose co mponent masses extend from 1.2 M & 
down to ~ 0.035 M flWhite et al. 1999D . these LMS and BD models give co- 
eval sequences for se veral other clusters like e.g. IC348, Taurus, Chamaeleon 
( Luhman et al. 20 03). Although the possibility that such an agreement be for- 
tuitous can not be totally excluded, the fact that these sequences are composed 
of hundred of sources lends confidence to the reliability of the models. 



3.2. Mass- magnitude relations 

Combining adaptative optics and accurate radial velocity surveys, Delfosse et 
al. (2000) and Segransan et al. (2003) have been able to determine mass- 
magnitude relationships (MMR) of nearby low-mass binaries down to the vicin- 
ity of the H-burning limit. Figure 01 displays the comparison of the Andersen 
(1991) and Segransan et al. (2003) data in the V band with different theoreti- 
cal MMRs, namely the parametrizations of Kroupa et al. (1993) (KTG), Reid 
et al. (2002) for My < 9 complemented by Delfosse et al. (2000) above this 
limit and the models of Baraffe et al. (1998) (BCAH98) for two isochrones. 
The KTG MMR gives an excellent parametrization of the data over the entire 
sample but fails to reproduce the flattening of the MMR near the low-mass 
end, whic h arises from the onset o f degeneracy near the bottom of the main 
sequence ( Chabrier Sz Baraffe 20001 ) . and thus yields a too steep slope. The 



Delfosse et al. (2000) parametrization, by construction, reproduces the data 
in the My=9-17 range. For My < 9, however, the parametrization of Reid et 
al. (2002) misses a few data, but more importantly does not yield the correct 
magnitude of the Sun for its age. The BCAH98 models give an excellent rep- 
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resentation for m ;> 0.4M©. Age effects due to stellar evolution start playing 
a role above m ~ 0.8 M©, where the bulk of the data is best reproduced for 
an age 1 Gyr, which is consistent with a stellar population belonging to the 
young disk (h < 100 pc). Below m ~ 0.4 M , the BCAH98 MMR clearly differs 
from the Delfosse et al. (2000) one. Since we know that the BCAH models 
overestimate the flux i n the V-band, due to still incomplete molecular opacities 
( Chabrier et al. 2000afi . mass function calculations (§4) have been conducted 
with the Delfosse et al. (2000) parametrization in this domain. The difference 
yields a maxim um ~ 16% discrepancy in the mass determination near My ~ 13 
( Chabrier 20021 . figure ^ displays a comparison of the data with the BCAH98 
calculations in all available observational bands. As seen in the figure, there is 
a remarkable agreement between the theoretical predictions and the data in all 
available infrared bands, the most favorable domain for LMS and BD observa- 
tions. Overall, the general agreement can be considered as excellent, and the 
error on mass functions derived when using these MMRs is smaller than the 
observational error bars in the LF. 




Figure 3. Comparison of various m-My relations 
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Figure 4. Compa rison of oberved i|Segransan et al. 200^ and theoretical 
l|Baraffe at al. 19981) mass-magnitude relationships in various bands 
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Figure 5. Comparison of the theoretical (Chabrier & Baraffe 2000) mass- 
radius relationship for low-mass stars with observations 



3.3. Mass-radius relation 

The radii of many LMS have been determined accurately from various tech- 
niques. Eclipsing binaries of course provide the most natural method but include 
only a limited number of systems below 1 Mq. Interferom etrv with the VLTI al- 
lows a precise determination of the radii of nearby binaries ( Segransan et al. 20 03). 



while transit observations from the OGLE microlensing survey improve signif- 
icantly the statistics. Figure [5] shows the comparison between the theoretical 
and observed radii from 0.8 Mq down to the HBMM, includi ng the determina - 
tion of the radius of the presently smallest H-burning object ()Pont et al. 2 005). 
The remarkable agreement between theory and observation gives confidence in 
the underlying physics used to determine the mechanical structure of these cool 
and dense objects, in particular the EOS (see §2). The observed radii for the 
eclips ing binaries, however, are found to be ~ 10% larger than the theoretical 
ones ([Torres fc Ribas 2002ft . As seen on the figure, however, all eclipsing bi- 
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nary radii but one lie above the other observed radii. It must be noted that all 
these eclipsing binaries are magnetically very active, a consequence of dynamo 
generated processes in objects at small orbital separation. It is thus natural 
to imagine that spot area cover a significant fraction of the irradiating surface 
of these objects, yielding smaller contraction along evolution. Although purely 
speculative at the present stage, this suggestion provides a viable explanation 
for the larger radius of eclipsing binaries compared to other binaries. 



4. Stellar and brown dwarf mass functions 
4.1. Galactic field 

A proper census of the number of stars and BDs bears important consequences 
not only for an accurate determination of the Galactic mass budget (about 70% 
of the Galactic mass is under the form of objects below a solar mass) but for 
our understanding of star formation. The determination of the stellar IMF can 
be done easily from observed luminosity functions (LF), using accurate mass- 
magnitude relationships (MMR). U sing the MMRs men tioned in §3.2, or fits 
based on observed nearby binaries (Del fosse et al. 200 0). the stellar IMF has 
been determined recently by Chabrier (2003, 2005), based on observed V-band 
and K-band LFs. Such an IMF, and comparison with other IMFs, is portrayed in 
Figure El The solid line superposed to the IMFs derived from the LFs illustrates 
the following analytical parametrization (in (log M ) _1 pc~ 3 ): 



r (log m — log 0.2) 2 i 
f(logm) = 0.093 xexp{- 1 % x (Q 55)2 }, ™ < 1 M 

= 0.041 m _1 - 35±0 - 3 ,m>lM Q (1) 

This IMF differs slig htly from the one derived in Chabrier (200 3), which 
was based on the 5-pc LF (|Dahn et al. 19861 : iHenry &: McCarthy 1990), whereas 
the present one is based on the revised 8-pc LF (lE,eid et al. 2002T ). The differ- 
ence at the low-mass end between the two parametrizations reflects the present 
uncertainty at the faint end of the disk LF, near the H-burning limit (spectral 
types ;> M5). Clearly a better determination of the faint part of the disk LF 
is required before the IMF can be determined with higher accuracy at the H- 
burning limit. Note that the field IMF is also representative of the bulge IMF 
(triangles), derived from the LF of Zoccali et al. (2000). 

The determination of the BD IMF is a complicated task. By definition BDs 
never reach thermal equilibrium and most of the BDs formed at the early stages 
of the Galaxy have now fainted to very faint luminosities, below present limits 
of detection. Observations are thus biased towards young and/or massive BDs. 
Monte-Carlo calculations, taking into account a mass probability distribution 
function P(m), given by the mass function, and a time probability distribution 
function, P(t), given by the star formation rate (SFR) must be conducted to 
obtain the BD IMF (assuming the mass- and age-distributions can be separated 
out). Such calculations have been carried out by Reid et al. (1999), Chabrier 
(2002, 2003, 2005), Burgasser (2004) and Allen et al. (2005). These calculations 
can be confronted to present observations of field BDs. 
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Figure 6. Disk IMF determined from various luminosity functions <f> 



Figure displays the calculated BD density distributions as a function of 
T e fj, L, Mk and Mj, obtained by aforementioned Monte Carlo calculations, us- 
ing the BD cooling models devel oped in the Lyon group, and the most recen t 
estimated LMS and BD densities (|Gizis et al. 200(1 : iBurgasser 20011 : ICruz 2 Q0j). 
These calculations include the correction for unresolved binaries. Binary frac- 
tions are based on the observed values i.e. a binary frequency decreasing from 
~ 50% to ~ 20% with decreasing mass from G-type stars to BDs. The agreement 
between the theoretical calculations and the observations is rather satisfactory, 
keeping in mind the remaining uncertainties both in BD cooling theory (in par- 
ticular in the late L early T transition domain) and in the observational determi- 
nations of the BD T e g, Mb i and number densities. It is interesting to note that 
the predicted dip around Mj ~ Mk ~ 13, M ho[ ~ 1 5 (IChabrier 2 003') has been 
confirmed by the most recent L-dwarf observations (|Cruz 20041 ). Thie dash-line 
displays the distributions obtained with a power-law IMF dn/dm oc m~ a with 
a = 1 , an upper limit for the coefficient ( Chabrier 20021 1 , with the same normal- 
ization at 0.1 Mq. As seen in the figure, only in the T dwarf domain do the two 
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Figure 7. Comparison of the observed LMS and BD distributions with the 
ones obtained with the IMF (1) (solid line) and with a power law IMF ^ oc 
m , with the same normalization at 0.1 Mq (dash- line). 



different IMFs yield significantly different densities, a point noted also by Allen 
et al. (2005). A more precise determination of the IMF thus necessitates large 
statistics in the T dwarf domain. This should be in reach with future dedicated 
large field surveys. A census of BDs, stars and stellar remnants in the Galaxy, 
based on these calculations, and their contribution to the Galactic mass budget 
is given in Chabrier (2003) with a slight revision in Chabrier (2005). Extend- 
ing IMF (1) down to one jupiter mass yields BD number and mass densities, 
ubd ~ n+/3 ~ 0.03 pc~ 3 and pbd ~ p*/30 ~ 1.0 x lO -3 M0pc -3 , respectively. 
The aforementioned power law IMF with a = 1 yields at least 5 times more 
BDs, depending on the minimum mass of the IMF. 

As shown in Chabrier (2003, Figure 10), the IMF (1) yields a mass-to- 
light ratio in various optical and infrared ba nds in very good agr eement with 
dynamical determinations in spiral galaxies ([Portinari et al. 20041 ). whereas a 
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Salpeter IMF extending down to the H-burning limit overestimates this ratio by 
a facto of 2. 

4.2. Young clusters 

As shown in Chabrier (2003, 2005) this field star+BD IMF (1) gives an excellent 
representation of the IMF of different young clusters with ages between ~ 1 Myr 
to ~ 100 Myrs, as derived from the observed LFs in various passbands, over the 
entire observed stellar and brown dwarf domain. Since these observations do 
not resolve multiple systems, the IMF in that case is the so-called system IMF, 
derived from the field IMF for resolved objects. The parametrization of this 
system IMF is given in Chabrier (2003, 2005). This similarity between field 
and cluster IMF clearly points towards a rather universal mechanism for star 
formation in various environments representative of the Galactic disk conditions. 
This issue is addressed in §6 

5. Present uncertainties in brown dwarf theory 

As outlined in the previous sections, the present LMS and BD theory has been 
succesfull in describing, or even in predicting the observed mechanical and ther- 
mal properties of these objects. Several shortcomings, however, remain and there 
is still room for improvement in various domains. The description of the spectral 
energy distribution of cool objects, in particular, needs further improvements. 

5.1. Young objects (~ 1-100 Myr) 

Comparison between observed and sythetic spectra shows an excess (lack) of flux 
in the J-band (H-band) for the latter ones, whereas phot ometry of optical color s 
shows that the models are too blue by ~ 0.3-0.5 mag (Cha brier et al. 20 00a). 
Both shortcomings point to the same kind of culprit, namely partially inaccurate 
molecular linelists. FeH lists, for example, cover only the 0.9-1.3 /im range. 
Water and TiO linelists and/or oscillator strengths are still imperfect as well. 
Not only these limitations affect the colors, in particular in the J and H bands, 
but they may also bear some impact on the atmospheric structures. Indeed, the 
derivation of fundamental parameters for young objects, i.e. objects still on the 
PMS, ranging from ~ 1 to ~ 600 Myr for low-mass stars, suggest some problems 
at young ages. Shortcomings in the theory are indeed more consequential for the 
evolution of still contracting objects, characterized by shorter Kelvin-Helmholtz 
time scales. Song et al. (2002) observed objects near the Lithium burning 
limit in a 12 Myr cluster and found no Lithium absorption in objects for which 
the models predicted a substantial amount of lithium. White & Hillenbrandt 
(2005) reach similar conclusions from their analysis of the optical spectrum of a 
spectroscopic binary with an accretion disk in Taurus Auriga. The age derived 
from lithium depletion is ~ 20 Myr whereas the age inferred from the isochrones 
at the distance of Taurus yields ~ 8 ± 3 Myr. The too low lithium-depletion 
in the models points to too hot T e g for a given mass and age (conversely the 
models suggest too small a mass for a given, observationnally-determined T e g). 

Caution must be taken, however, before drawing robust conclusions from 
such an observational analysis. Indeed, the Song et al. (2002) analysis is based 
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on T e ff-color relations characteristic of main sequence objects (|Bessel 199lh and 
uncertainties reach at least ±100 K. As for St 34, the object of White & Hillen- 
brandt (2005), a subse quent spectroscopic analysis with the Spitzer Telescope 
((Hartmann et al. 20051 ) shows that the object is more likely an older object, 
possibly as old as 25 Myr, in the foreground of Taurus. The ages derived from 
the models using the lithium depletion or the isochrones would then be fully 
consistent. Furthermore, all these young objects show H a emission lines, a sign 
of accretion and/or activity. Intense H Q emission implies the presence of a 
hot (> 10 4 K), ionizing flux which may penetrate at least the outer parts of 
the photosphere, where the cores of Li I resonance lines are formed. This will 
yield additional (UV) continuum and overionization and thus reduced equivalent 
widths of absorption lines. Indeed, Zapatero-Osorio et al. (2002) and Kenyon 
et al. (2005) found significant dispersion in the lithium equivalent widths of 
M stars for spectral types cooler than M3.5 in the Sigma Ori cluster. Lithium 
abundance analysis in the presence of strong H a emission must thus be consid- 
ered very cautiously, as what could be considered as l ithium depletion mig ht in 
fact be due to non-photospheric line formation effect (|Pavlenko et al. 19951 ). 

In a different analysis, Mohanty et al. (2004) suggest also some incon- 
sistency in the theoretical mass-radius relationship between the determination 
based on spectral analysis and the one obtained from evolutionary models. Al- 
though, as pointed out by these authors themselves, gravity and age determina- 
tions from the observed spectra remain rather uncertain for cool objects where 
the atmosphere is dominated by molecular lines, their careful analysis suggests 
too small a gravity (i.e. mass) for a given T e g at a given age, in particular for low- 
gravity objects. The same conclusion is obtained from the recent observations 
of AB Dor C, a low-mass companion to the young (~ 30-100 Myr) star AB Dor 
A, by Close et al. (2005), namely too hot T e g, and thus too bright luminosities 
from the models for a given mass and age at young ages. The T e g determination 
of Close et al. (2005), however, is based on an M dwarf T e fj-Spectral type (Sp) 
calibration which is not valid for younger, lower gravity objects like AB Dor C . 
Moreover, anew analysis of the age of the AB Dor system d ammar, et al. 2005a l ) 
shows that the observations are indeed in fairly good (la) agreement with the 
model predictions. The same conclusion is reached if AB Dor C is an unresolved 
BD binary system itself, in which case the obse rved T P ff and magni tudes would 
superbly corroborate the theoretical predictions ( Marois et al. 20051 ). Radial ve- 
locity analysis or high resolution images should eventually settle this question. 

This illustrates how much care must be taken from observational analysis 
before reaching robust conclusions about the validity of theoretical models. T e g 
determinations from spectral types for young (low gravity) objects, in particular, 
remain very uncertain. Indeed, as shown by Luhman (1999), young LMS and 
BD have spectral types which lie between the M dwarf scale and the giant scale 
and so far no robust T e g-Sp relation exists for objects with in-between gravi- 
ties. Clearly, the determination of fundamental parameters from observations of 
young objects remains hampered by substantial uncertainties and, although the- 
oretical models at young ages certainly need to be improved, conclusions from 
observational analysis must be considered with due caution. 
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5.2. Very young objects (<; 1 Myr) 

Baraffe et al. (2002) have done an extensive analysis of uncertainties in the mod- 
els for very young ages 1 Myr). As shown by these authors, unknown initial 
conditions yield different evolutionary paths until at least ~ 1 Myr and mass-age 
determinations from comparison between observations and evolutionary tracks 
is totally meaningless in this case. Other processes, like e.g. the efficiency of con- 
vection both in the interior and the atmosphere of low gravity objects remains 
totaly undetermined and can modify appreciably the PMS contraction sequence. 
Evolution, for example, does not necessarily proceed at constant T e g for a given 
mass, as is th e case only for fully adiabatic convection in presence of di ssoci- 
ation regions (Havas hi 196 it IHavashi k, Nakano 19631 : iBaraffe et al. 20 02). For 



this reason, our group decided not to publish models for such very young ages. 
It is worrying, however, that some observers use blindly available models at such 
ages although at best these models are as inaccurate as ours (see e.g. Fig. 8 of 
Baraffe et al. (2002)), and in some cases are based on simplified atmospheric 
structures and remain inaccurate at any age ! Unfortunately, determinations 
of young cluster mass functions based on such erroneous analysis are claimed 
periodically. We can not stress enough the necessity to consider such determi- 
nations not seriously and to remain aware that many claimed low-mass IMF 
determinations for young clusters in the litterature have no reliability ! 



6. Brown dwarfs versus planets. Brown dwarf and star formation 

The distinction between BD and giant planets has become these days a topic of 
intense debate. In 2003, the IAU has adopted the deuterium-burning minimum 
mass, rriDBMM — 0.012 M & (Saumon et al. 1996, Chabrier et al. 2000b) as 
the official distinction bet ween the two types of objects. We ha ve discussed this 
limit in previous reviews (IChabrier 20031 : Ichabrier et al. 2004ft and shown that 



it does not rely on any robust physical ground and is a pure semantic definition. 
Such an approach is certainly not satisfactory from the scientific point of view 
and should be abandoned, for it brings a lot of confusion in the field. This dis- 
tinction is based on a model of star formation where deuterium-burning plays a 
central role (Shu et al. 1977). Indeed deuterium burning generates a convective 
instability which in turn leads to the generation of a magnetic wind through dy- 
namo generation. This magneto-centrifugal wind is necessary to sweep away the 
surrounding accreting gas and thus leads to the formation of a star/BD embryo. 
This scenario can now be abandoned for various reasons. First of all observa- 
tions of star forming regions show that star formation occurs in general on a time 
scale much shorter than the standard ambipolar diffusion time scale character- 
istic of the previous scenario, namely a few dynamical or turbulent time scales 
(see e.g. Hartmann et al. (2001), Kirk et al. (2005) their Fig. 11). Second, 
all BD interiors are convective, with or without deuterium burning. Although, 
admitedly, this holds for a non-accreting, or at least non spherically accret- 
ing protostar and has not been demonstrated to remain true for an accreting 
object. The phase of spherical accretion, the so-called "c lass 0" phase, how - 



ever, is observed to take place over a very short time scale (| Andre et al. 2 000). 
Very rapidly accretion is taking place throug h a disk or along magnetic poles 
( Bontemps et al. 19961 : iHenriksen et al. 1997t ) and thus the internal structure 
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of the protostar very likely remains convective ((Hartmann et al. 1997h . Finally, 
even though magnetic winds may play a role in star/BD formation, it is unlikely 
that they play a dominant role and gas reservoirs limited by turbulence are ap- 
pealing alternative mechanisms to limit the accretion on the central core. In 
any event, as mentioned above, deuterium burning can not be considered as a 
necessary mechanism to trigger star formation and thus can not be considered 
as a robust limit for star formation. In fact the observation of free floating ob- 
jects wit h masses of the or der of a few jupiter masses in (low extinction) young 
clusters (Bciar ct al. 200l]) shows that star and BD formation extends down to 
jupiter-like masses. 

The ~10 jupiter-mass limit, however, illustrates also the so-cal l ed opacity - 
limited minimum mass for fragmentation ( Low Sz Lvnden-Bell 19761 : ISilk 19771 ) . 
This limit is generally considered as robust because of the weak dependence of 
this minimum mass upon (dust) opaci ty, m m j n oc k 1 / 7 . Effects like magnetic 
field dBoss 200lft . anisotropy or shocks (|Bovd Whitworth 20051 ). however, can 
significantly lower this limit, down to a few jupiter masses. 

We thus suggest the following clarification for "star", "brown dwarf, "sub 
brown dwarf, "planetary object" or "planet" denominations. Objects formed 
either in systems or as isolated objects from the collapse of a cloud are either 
stars or brown dwarfs. The physically grounded limit between the two types of 
object s is the hydrogen-burning minimum mass, m ~ 0.07 M@ for solar abun- 
dance ((Chabrier &: Baraffe 19971 ) . below which the object can not reach thermal 
equilibrium and keeps contracting steadily with time (Chabrier & Baraffe 2000 
Fig. 2). Given this, there are deuterium-burning brown dwarfs and non deu- 
terium burning brown dwarfs, depending whether their mass exceeds or not the 
deuterium-burning minimum mass, exactly as stars above ~ 1.5 Mq ignite the 
CNO burning cycle in their core and stars below this mass do not. Both types 
of objects, however, are called stars and not "substars" or other exotic name ! 
Planets, on the other hand, are objects formed in a protoplanetary disk around 
a parent star, with a high mass ratio. They should thus have enhanced average 
abundances of heavy elements. Although the direct spectroscopic observational 
confirmation of such a diagnostic remains for now out of reach, the recent ob- 
servation of the transiting planet HD 149026B, with an inferred 70 of heavy 
element for a planet mass m sini = 0.36 uij up (Sato et al. 2005) seems to com- 
fort this suggestion. Therefore, the various "direct" detections of exoplanets at 
large orbital distances (^ 50 AU) claimed recently in the literature, and based 
on the aforementioned IAU definition, are most likely in reality wide binary 
brown dwarfs (see Gizis et al. (2005)), by itself an interesting result. 

This brings us to the issue of star and BD formation. As mentioned above, 
the conventional scenario for star formation, which assumes that proto-stellar 
cores are in quasi-static equilibrium against gravitation and form stars once 
magnetic support is lost through ambipolar diffusion ( Shu et al. 1977T ) . can 
now be abandoned on various grounds. The modern picture of star forma- 
tion suggests that cores a re formed by compressible motions in the turbulent 
velocity field of the cloud ( Ballesteros-Paredes et al. 20031 : iKlessen et al. 2005; 
iPadoan Sz Nordlund 2002T ). Indeed, large scale (> 0.1 pc) cloud structures are 
dominated by (supersonic) turbulence. At small (stellar) scales, gravity takes 
over and fragmentation becomes gravitational rather than turbulent, with only 
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Figure 8. Comparison between the system IMF representative of the Galac- 
tic disk and young clusters (Chabrier 2003) and the distribution of core masses 
obtained from compressible MHD turbulence (Padoan & Nordlund 2002) 



the densest cores with gravitational energy exceeding the local internal energy 
collapsing into stars. In this scenario, progressively smaller density peaks con- 
tain progressively smaller mass fractions, yielding de creasing star formation effi- 
ciency with decreasing scales b elow the Jeans mass (|Padoan fc Nordlund 20021 : 
IVazauez-Semadeni et al. 2003t h In these cal culations, BDs represent si mply the 
low-mass part of the star formation process (|Padoan fc Nordlund 2004 ). Figure 
|S] compares the mass spec trum obtained from MHD simulations of supersonic, 
super Alfvenic turbulence ( Padoan Sz Nordlund 200 2) with the system IMF ob- 
tained for the Galactic disk and young clusters in §4 (indeed, the identified 
prestellar cores in the simulation correspond to multiple systems, for some of 
these cores might eventually become unstable and fragment into more individual 
objects). Although such a comparison must be considered with caution before 
drawing any conclusion (e.g. gravity, which will be dominant at the stellar scale, 
is not included in the simulations yet), the at least qualitative and even semi- 
quantitative agreement between the simulations and the IMF representative of 
the field and young clusters is encouraging. This shows that BDs form in adequat 
numbers from the same fragmentation mechanism as for star formation. Various 
observations of disk accreting BDs show indeed that BDs and stars form from the 
same underlying mechanism fwhite &z Basri 20031: iJavawardhana et af~ 2003; 
iLuhman et al. 2005bl : iGizis et al. 20051 : iKenvon et al. 20051 ). Note that in these 
simulations of supersonic turbulence, only a few percents of the total mass end 
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up into the collapsing cores after one dynamical time, providing a solution to the 
problematic high efficiency problem associated with turbulence-driven star for- 
mation. Although admitedly time scale may remain an issue. This new picture 
thus combines turbulence, as the initial driving mechanism for fragmentation, 
and gravity, providing a natural explanation for a (scale free) power-law IMF 
above a critical mass, namely the mean thermal Jeans mass (mj), and a lognor- 
mal distribution below, due to the fact that only the densest cores, exceeding 

the local Jeans mass, will collapse into bound objects. 

The a l ternative scenario of formation of BDs by ejection ( Reipurth Sz Clarke 200ll : 
iBoss 20011 : iBate et al. 2 002) can now be abandoned as a dominant scenario for 
BD formation for various reasons. First of all, the aforementioned observa- 
tions of disk accretion in BDs exhibit characteristics similar to the more mas- 
sive T Tauri st ars. Second of all, the recent observation s of wide binary BDs 
( Luhman 20041 : iForveille et al. 20041 : IChauvin et al. 20051 : iBillieres et al 2005) 
contradict the predictions of such a scenario as a dominant formation mech- 
anism for BDs. At last, the observed high binary frequency of VLMS/BDs, 
in particular for close binaries, rules out at >99 confidence level a dynamical 
origin for VLMS/BDs, as obtained from N-body models (which neglect in partic- 
ular gas-dynamic harden ing mechanisms) or hydro dynamical SPH simulations 
( Maxted fc Jeffries 2005h . 



7. Conclusion and perspectives 

As mentioned in the introduction, accurate models for low-mass stars, brown 
dwarfs and giant planets are needed to shed light on the observable properties 
of these objects and to provide guidance to the ongoing and future surveys aimed 
at understanding their formation, structure and evolution and at revealing their 
contribution to the Galactic population. Tremendous progress has been realyzed 
from this point of view within the recent years. As outlined in this review, the 
present modern theory of LMS/BD structure and evolution has been very suc- 
cesfull in accurately predicting several stringent observational constraints : spec- 
tra, mass-magnitude and mass-radius relationships, color-magnitude diagrams. 
Any theory aimed at describing the mechanical and thermal properties of these 
objects must be confronted to these experimental/observational constraints in 
order to assess its degree of validity. Improvement in the theory is still neces- 
sary, in particular for young, low gravity objects, and will proceed along with 
the discovery of many more substellar objects, down to Jupiter- like masses. 

Present determinations of the star+BD IMF give satisfactory comparisons 
with observational determinations and yield reasonably robust number-density 
and mass-density estimates for the stellar and BD census in the Galaxy. The 
increasing number of observed BDs will improve the accuracy of these deter- 
minations. Although such an IMF can be described by a series of power laws, 
such an approach, though convenient, is not satisfactory form the physics point 
of view as it can not be rooted to any physical mechanism for star formation. 
Indeed any change in the slope implies a different characteristic mass. A gen- 
eral Salpeter-like power law above the average thermal Jeans mass (~ IMq for 
present day conditions) rolling down into a lognormal form below this mass, on 
the contrary, is consistent with physical mechanisms like turbulence-driven frag- 
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mentation. Although a complete theory for star/BD formation is still elusive, 
turbulence fragmentation at large scale indeed provides an appealing solution for 
the "universal" triggering mechanism and leads to the formation of stars and 
BDs from the same initial process. Processes like magnetic winds, accretion, 
dynamical interactions will certainly play some role in shaping up the final IMF 
but they do not seem to play a dominant role in the underlying general process. 

At last, we argue that the deuterium-burning limit as a definition of BD 
vs planet distinction should be abandoned, for it does not fit into a modern 
approach of star formation, and BD formation very likely extends below this 
limit. 
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